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Abstract 

It is known that after Affleck-Dine baryogenesis, spatial inhomogeneities of Affleck- 
Dine field grow into non-topological solitons called Q-balls. In gauge mediated SUSY 
breaking models, sufficiently large Q-balls with baryon charge are stable while Q-balls 
with lepton charge can always decay into leptons. For a Q-ball that carries nonzero 
B and L charges, the difference between the baryonic component and the leptonic 
component in decay rate may induce nonzero electric charge on the Q-ball. This 
implies that charged Q-ball, also called gauged Q-ball, may emerge in our universe. 
In this paper, we investigate two complex scalar fields, a baryonic scalar field and a 
leptonic one, in an Abelian gauge theory. We find stable solutions of gauged Q-balls for 
different baryon and lepton charges. Those solutions shows that a Coulomb potential 
arises and the Q-ball becomes electrically charged as expected. It is energetically 
favored that some amount of leptonic component decays, but there is an upper bound 
on its amount due to the Coulomb force. The baryonic decay also becomes possible 
by virtue of electrical repulsion and we find the condition to suppress it so that the 
charged Q-balls can survive in the universe. 



1 Introduction 


Affleck-Dine mechanism [1] is a promising candidate for baryogenesis in supersymmetric 
(SUSY) theories due to its consistency with the observational bound on reheating tempera¬ 
ture which avoids the gravitino problem [2], In the Affleck-Dine mechanism, baryon asymme¬ 
try is generated through the dynamics in the phase direction of a complex scalar field which 
carries nonzero baryon charge [3]. The scalar field is called Affleck-Dine field. Affleck-Dine 
field is spatially homogeneous when it starts to oscillate, but spatial inhomogeneities due 
to quantum fluctuations grows exponentially into non-topological solitons, which are called 
Q-balls [4, 5, 6]. Q-ball is a spherical condensate of a scalar field and is defined as a solution 
in a global U(l) theory which minimizes energy of the system with its charge fixed [7]. Q- 
ball is known to decay into quarks or leptons, so that the final baryon number is carried by 
quarks produced through the decay of Q-balls. However, in gauge mediated SUSY breaking 
models, a baryonic Q-ball, with sufficiently large charge, can be stable against decay into 
nuclei when energy per charge of Q-ball is smaller than the proton mass [8]. On the other 
hand, for a leptonic Q-ball, there exist decay channels into leptons. In this case, baryon 
number of the universe is generated from lepton asymmetry through sphaleron effect, i.e. 
leptogenesis. 

We focus on Q-balls that carry both baryon and lepton charges. In fact, these Q-balls can 
be formed when we consider the Affleck-Dine baryogenesis with u c u c d c e c flat direction, for 
instance. In this case, it is possible that its lepton component can decay into leptons while its 
baryon component cannot decay into baryons. This implies that the difference between the 
baryonic component and the leptonic component in decay rate may induce an electric charge. 
Therefore, through the decay of the leptonic components only, an electric charge is induced 
even if the neutral Q-ball was formed at the beginning. The electric charge is expected to 
make differences in the experimental signatures of the relic Q-balls. For instance, the neutral 
Q-balls can be detected by Super-Kamiokande [9, 10, 11], which probes the absorption of 
protons, but this detector is not suited for detection of the charged Q-balls since the charged 
Q-balls cannot absorb the protons due to the electrical repulsion. The charged Q-balls are 
likely to behave as some kind of nuclei, and are known to be detectable by such detectors as 
MACRO [9, 12] and the observational bounds on mass and flux of the relic charged Q-balls 
are obtained [10]. 

In fact, electrically charged Q-ball has been studied in the literature and is called gauged 
Q-ball [13] . However, the previous works studied gauged Q-balls in one scalar held theories, 
so that their results cannot be applied to Q-balls generated after the Affleck-Dine baryoge¬ 
nesis. There is a previous work which also discussed the evolution of Q-balls in two scalar 
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model [14, 15], in which the gauge field was neglected. However, in the course of decay, it 
is expected that an electrical repulsion arises, so that the gauge field must be taken into 
account. 

In this paper, we consider a simplified model where there are two complex scalar fields 
and U(l) gauge field. One of the complex fields carries baryon charge and positive electric 
charge, while the other one carries lepton charge and negative electric charge. Our main 
purpose is to demonstrate that gauged Q-balls may be realized in our universe through the 
decay process of the lepton component, even if Q-balls are initially neutral. We find stable 
solutions for different baryon and lepton charges, taking into account the effect of the gauge 
field. If we suppose that only the lcptonic component decays off, a sequence of solutions 
with B = const, represents the decay process. We examine quantitatively whether such a 
process is energetically allowed. 

In the following sections, we first review the main properties of global Q-balls and gauged 
Q-balls, i.e. neutral and charged Q-balls respectively. Subsequently, in order to discuss the 
evolution of a Q-ball which is formed from the flat direction, we consider a two scalar model 
and find sequences of gauged Q-ball solutions with B = const. We calculate the total energy 
of Q-balls and gauge fields and examine whether the leptonic decay is energetically favorable. 
Finally, we give our conclusions in Sec. 5. 


2 Global Q-ball 


In this section, we briefly review the main properties of a global Q-ball, which is a stable 
configuration of a complex scalar field with a fixed conserved charge. Consider a theory of 
a complex scalar field with a global U(l) charge. The Lagrangian is written as 


C = - F(4>) 


( 1 ) 


where V($) is a scalar potential and is normalized so that V(0) = 0. Let us renormalize the 
field for later convenience as 



( 2 ) 


The global U(l) charge density q is given by 



( 3 ) 
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Q-ball is defined as a solution which minimizes energy of the system with its charge fixed. 
Using the Lagrange multiplier method, we only need to minimize the following function: 


where E is given by 


E u = E + to 



E= (Ex 


^ (I0I 2 + |V0H + 


Equation (4) is rewritten as 


Eu = J d 3 x^\d t tj> - | 2 + 

VM) = V(® - ^ 2 |«f. 


(TX 


i| V4,\ 2 + VM 


+ uQ , 


(4) 

(5) 

( 6 ) 
(7) 


By minimizing the first term of Eq. (6), we can derive the time dependence of the solution 
as 


(t>(x, t) = e lu}t (p(x). 


( 8 ) 


Moreover, it is known that the solution which minimizes the second term of Eq. (6) is real and 
spatially symmetric [16]. Therefore, the radial direction 0(r) is a solution of the following 
equation: 


d 2 (j) 2 d(f> 
dr 2 r d r 


cu 2 (f)(r) 


dV(</>) 

dcf) 


0. 


(9) 


Here, in order to avoid a singularity at r = 0 and to find a specially localized solution, we 
set boundary conditions as 

^(0) = 0, <t>( oo) = 0. (10) 

Now let us find the condition that there exists a solution of Eq. (9). First, we redefine 
the spatial coordinate and the potential as 


r —>■ t, 

-VM = \uj 2 ^ - 


The equation of motion then becomes 


+ - 


2 ; dV,„ 


t dcj) 


= 0 


(11) 

( 12 ) 


(13) 
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Figure 1: Several kinds of shape of —V u ((p) depending on c o. (I) : cUq < u 2 < E"(0), (II) : 
co 2 < cUq, (HI) : co 2 > E"(0). The potential of type (II) is appropriate for the existence of 
the solution. 


and we see that it is equivalent to the equation of motion of a classical point particle under 
the potential —V u . The boundary conditions Eq. (10) imply that the classical particle is 
initially at rest and converges toward the origin asymptotically. In order that this kind of 
motion is possible, the potential — should be in an appropriate shape as shown in Fig. 1. 
Quantitatively, — must satisfy the following conditions: 

ma x.[-VJ<j>)] > 0, (14) 


- 00 ) < 0 . 


(15) 


The first condition is necessary in order for the particle to possess the potential energy 
enough to reach the origin. The second means that the particle must be subjected to a 
backward force near the origin and is also necessary in order for the particle to stop at the 
origin, not to roll over it. The conditions above are rewritten as the conditions to u as 


(Ug < oj 2 < y"(o), 


Wo = min 


0 2 


(16) 

(17) 
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Furthermore, it generally holds that 


dE 

^ =UJ i 


dQ 


(18) 


for any Q-ball solution, which can be easily shown by taking a variation of energy in the 
following way: 



(19) 


where we used the equation of motion Eq. (9). Therefore, the condition Eq. (16) is rewritten 
as 


dE /— 

Uo< dQ < = m». 


( 20 ) 


It is noted that the second inequality indicates that a Q-ball with a charge Q is energetically 
favorable compared to a Q-ball with a charge Q — 1 and a particle, which is consistent with 
the definition of Q-ball solution. 

Here, let us identify the scalar held 0 as a D- and F- hat direction in gauge mediated 
SUSY breaking models. In this case, the potential of the hat direction is mainly given by its 
soft mass term. However, the soft mass is suppressed for energy scale larger than messenger 
scale and the potential becomes hat 1 . If we approximate the potential as V — Vq — const., 
there exists an analytic solution [8] given by 



( 21 ) 


( 22 ) 


whose energy is then written as 



(23) 


Then, from dE/dQ = oj, which is proven above, 


0 1 = \/2nV,l' 4 Q- 1/a . 


(24) 


1 The exact form of the potential is derived in Ref. [17]. 
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This indicates that for a large charge, dE/dQ may become small enough. Therefore, for a 
baryonic Q-ball with a large baryon number, it may hold that 


dE 

dB 


< rn p , 


(25) 


where m p (~ lGeV) is the proton mass. This means that the Q-ball is stable against the 
decay into protons. On the other hand, if Q-ball carries a lepton number, it can decay into 
leptons. This is the motivation of our assumption in Sec. 4 that baryonic component of 
Q-ball is stable while leptonic component can decay into leptons. 


3 Gauged Q-ball 


In this section, we consider a complex scalar field that is charged under U(l) gauge symmetry. 
Although this toy model is not motivated in SUSY theories, we investigate it to see the effect 
of gauge force on Q-balls. In this case, we must solve the equation of motion for the gauge 
field A p as well as the complex scalar field. The spatially localized configuration in this 
theory is called gauged Q-ball. 

The Lagrangian is written as 

C = - U(4>) - (26) 

D p = d p - ieAp (27) 

where V ($) is a scalar potential. We parameterize the scalar field in the same way as in the 
case of global Q-ball: 

(28) 

4>(x,t) = e lut 4>(r ). (29) 

For the gauge field, we adopt the following parameterization [13]. 

4) = A 0 (r), (30) 

Ai = 0. (31) 


The first indicates that we are searching the spatially symmetric solution and the second 
implies that we are assuming the absence of a magnetic field, which means in turn the 
absence of an electric current. The equations of motion are then given by 


d 2 <j) 
dr 2 

ii 

dr 2 


2 d<p 
r dr 
2dg_ 
r dr 


+ 0(] z - — 


dV 

dcf) 



0, 


(32) 

(33) 
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where we redefine the gauge field to absorb oo as g = oo — eA 0 . Note that g is gauge invariant. 
We set boundary conditions as 


0(°°) — 0, ^(0) = 0, (34) 

g( oo)=w, |(0) = 0, (35) 

to avoid singularities at r = 0. Note that when 0(r) —* 0 as r —» oo, the gauge field 
asymptotes to a certain constant as r —» cx) by Eq. (33). Thus the boundary condition 
g( oo) = a; is just a definition of a;. Note that Eq. (33) can be rewritten as 

(■ r 2 g')' = e 2 r 2 (j) 2 g. (36) 


This implies that if g( 0) > o, then g' becomes positive for r > 0, so that g increases, while 
in the opposite case g(0) < 0, g decreases. In either case, g 2 always increases. 

Here, let us consider an analogy in a similar way as considered in the previous section 
(Fig. 2). Eq. (32) is analogous to the equation of motion for a classical particle which is 
subjected to the potential —V g = — V + \g 2 <\> 2 ■ As mentioned above, g 2 always increases, 
so that the effective mass term g 2 ^ 2 /2 increases with time. From this we can derive some 
important properties of gauged Q-ball. Since the effective mass increases, it is energetically 
possible that the particle reaches the origin even if it moves away from the origin at the 
beginning, which means that there exists radially non-monotonic solution. We can interpret 
this kind of solution as the result of the scalar field being pushed outward due to the electrical 
repulsion. We show both kinds of solutions for gauge mediation-like model in Fig. 3, where 
we approximate the potential as H( < h) = m|ln(l + <& 2 /m\). Indeed the non-monotonic 
solutions arise for charges larger than those of the monotonic ones. 

It holds that 



(37) 


for gauged Q-ball as well, whose proof is similar to that for global Q-ball [18]. The energy 
of the Q-ball is given by 


E 


(TX 


-(V</>) 2 + -(VA 0 ) 2 + -0 2 (w — eA 0 ) 2 + V(4>) 


(38) 
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Figure 2: Mechanical analogy of gauged Q-ball 
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Figure 3: Gauged Q-balls in gauge mediation-like model with e 2 = 0.002 and = 1. 
whose variation leads to 


5E 


J d 3 x [—Acj)8(f) — Aq5/S.Aq + (j)8(j)(u! — eA 0 ) 2 + (u 

j d 3 x [2 cf)8(f>(uj — eA 0 ) 2 + (u> — eA 0 )8(uj — eA 0 )(p 2 

J d 3 x [(ca — eA 0 )8q — AohAA 0 ] 

cu8Q — J d 3 x [eA 0 8q + A 0 5AA 0 \ 
u5Q. 


eA 0 )8(ui — eA 0 )(j) 2 + V'5<ft\ 
AohAkLo] 


( 39 ) 














log(e 2 ) 


Figure 4: Illustration of upper limits on Q, and where dE/dQ = m$, 


Here we used Eqs. (32) and (33), and the charge of the Q-ball is given by 

Q =h d3x 

= / d 3 x(w — eA 0 )(j) 2 . 


(40) 

(41) 


The charge of a gauged Q-ball has an upper limit, above which there is no localized 
solution. This is contrast to the case of global Q-balls, where Q-ball solutions exist for 
arbitrarily large Q. The blue circles in Fig. 4 are the upper limits on charges, and the results 
can be fitted as 


logio <?«!»* = -1.0 x 10 1 - 1.9 log 10 (e 2 ). 


(42) 


which is shown by the blue line. 

For a global Q-ball, the condition 


dE 

— < m$, 


(43) 


was necessary for the existence of solution (see Eq. (16)). However, for a gauged Q-ball, it 
is possible that 


dE 

dQ 


> m$, 


(44) 


9 




which implies a Q-ball with a charge Q — 1 and a particle at infinity are energetically favorable 
compared to a Q-ball with a charge Q. This behavior arises when the charge is large enough 
and therefore may be interpreted as a result of electrical repulsion. Even so, a gauged Q-ball 
may exist as itself because a Q-ball with a charge 0 is still energetically favorable compared 
to a Q-ball with a charge Q — 1 and a particle near the surface just after emission, due to 
the large Coulomb potential. Thus, gauged Q-ball is expected to be a metastable solution if 
Eq. (44) is satisfied. The red circles in Fig. 4 demonstrates when dE/dQ becomes m$. The 
results can be fitted as 


logioQ = -11 x 10 1 - 1.61og 10 (e 2 ). 


(45) 


This dependence can be explained in the following way. If we approximate the energy of an 
emitted particle dE/dQ as the energy with electricity switched off, plus Coulomb energy, 


dE 

dQ 


w 0 + 


Wo + 


Wo + 




47r R 


47r 2 


e 2 Q 


w 0 


V 4 n—I / 4 


2 \/ 2 ' 


7T 


r„'*q 


(46) 


where we used the analytic expression in the previous section with R = n/uJo, ujq = 
V2-kVo II Q^ 1/a . The charge of Q-ball at which dE/dQ = m<s> is thus given by 


Q = 


f 2y/2n m$\ 


(47) 


or 


logio <2=g logic 


/ 2y / 2vrm$ \ 4 2 

- o lo gio(e ) 


\ vE ) 


~ 9.9 x 10 


-l 


1.3 log 


3 

,2\ 


101 


(48) 


where we used coo and set Vo into = 1.6m$, whose value is appropriate for the 

solutions we are dealing with. This estimation roughly explains our numerical solutions of 
Fig. 4. 


4 Q-balls in two scalar model with U(l) gauge field 

Since our main interest is in the evolution of global Q-balls formed from the flat direction 
which is responsible for the Affleck-Dine mechanism, we must consider the case of several 
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numbers of scalar fields coupled with a gauge field. Here we consider the simplest case in 
which the flat direction consists of two scalar fields which carry baryon and lepton numbers 
respectively, and the gauge held is Abelian. The electric charges must be opposite since the 
hat direction is neutral. The Lagrangian is then written as 


£ = (D^YD^ + (D,$ 2 )*D^ 2 - $ 2 ) - 

D^ 1 = (d tl - ieAp)®!, 

D^ 2 = (d/j, + ieA^)<& 2 

and baryon and lepton charges are 


1 

i j 

[ = 

j d 3 xb , 

1 

* J 

f d*x(%D^2 ~ $2(A$2)*) = j 

j d 3 xl 


(49) 

(50) 

(51) 


(52) 

(53) 


where b and l are baryon and lepton number densities. Since we assign the positive and 
negative electric charges for B and L components, respectively, the total electric charge is 
given by 


Q = B-L. 


(54) 


We find stable solutions and calculate their energies, and examine if leptonic decay is 
energetically allowed. First, we adopt the same parameterization as before: 


1 = !> 2 

(55) 

M X A) = e twit (j) i(r), 

(56) 

<MM) = e^fair), 

(57) 

Ai = 0, 

(58) 

bio — A)( r )- 

(59) 


The equations of motion then become 


d 2 <P i 2 #i , ( , , 2 d V_ 

H-j-b — eAo) — —— — 0, 

r dr ocpi 


dr 2 

d 2 (p 2 2d(t> 2 2 

, n H-J-h 0 2(^2 + CAq) 

dr z r dr 

d 2 A 0 2dA 0 2 / 

+ ^ (W1 


eA n 


dV 

tt-=0, 

d(j) 2 

- e<p 2 (uj 2 + cAq) 


= 0 


(60) 

(61) 

(62) 
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with the boundary conditions given by 


We prove here 


0i(oo) = 0 2 (oo) = o, 

^(0) = ^(0)=0, 

dr dr 

4)(oo) = ^h(O) = 0. 

dr 



(63) 

(64) 

(65) 

( 66 ) 

(67) 


for later use, which is analogous to dE/dQ = u for 1-scalar gauged Q-ball. The energy of 
the system is 


E = d 3 x 




e^o) 2 + + e^o) 2 + fd(0 1; 02) 

( 68 ) 


and its variation with respect to 0i, 0 2 and A 0 is given by 


/OT T OT T 

d 3 X — A01<501 — A(j) 2 5(J) 2 — AM^O + _ 501 + _ (502 

<90 1 C/02 

+ J d 3 x [0i<50i((Ui — eA 0 ) 2 + 0i(ci/i — eAo)(f(ct/i — eAg)] 

+ J d 3 X [02(502(^2 + cAq) 2 + 02(o/ 2 + e -A())5(ti/2 + C-A)] 

= j d 3 x [(o/i — eA 0 )5b + ( 0/2 + eA 0 )5l — A 0 5AA 0 ] 

= oji5B + oj 2 SL + J d 3 xA 0 [e(—5b + Sl) — SAA 0 ] 

= ojiSB + oj 2 5L, 


(69) 


where we used Eqs. (60), (61), and (62). The proof can be easily generalized into the case 
of arbitrary number of scalar fields. 

We assume that the decay of the Q-ball always takes place if energetically allowed and 
that the evolution of the Q-ball can be approximated as a sequence of gauged Q-ball solutions. 
Then the decay of the leptonic component only can be represented by a sequence of gauged Q- 
ball solutions with the same baryon number arranged in descending order of lepton number, 
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which is expected to decrease due to the decay. The results for B = 1.7x 10 4 and B = 8.4x 10 4 
are shown in Fig. 5 and in Fig. 6 respectively, where we also used the following approximate 
form of the gauge mediation potential: 

ln (l+^) + m% 2 In ^1 + ^ + y (|$i| 2 - |$ 2 | 2 ) 2 ■ (70) 

Here, we include the D-term potential, which arises since the D-flat condition |$i| = |<T> 2 1 
is not valid anymore, and in addition, we assume m^ 1 = m$ 2 = m$. We see that as the 
leptonic component decays, the gauge field, or the Coulomb potential arises. Therefore, it 
is expected that through decay process, the initially formed neutral Q-ball may evolve into 
the charged or gauged Q-ball, which means that the charged Q-balls may emerge in our 
universe. Next, if we look at the energy, we see that in Fig. 5, the energy decreases along 
the decay, which indicates that the particle which comes out from the Q-ball has positive 
energy. This means that the Q-ball emits free (i.e. unbounded to the Q-ball) particles until 
the leptonic component completely vanishes. Whereas in Fig. 6, the energy starts to increase 
in the middle of the decay, which means that the energy of the emitted particle becomes 
negative. This in turn means that the particle starts to be bound to the Q-ball. This may be 
understood in the context of quantum mechanics of many-body system as a cloudy bound 
state which is analogous to that of an atomic system. We also see a slight decrease in energy 
near the end of the decay, which means that the emitted particle becomes free again. This is 
interpreted as follows. As shown in Fig. 6, the leptonic component concentrates at the center 
while the baryonic one locates outside away from the center. Since the baryonic component 
is far enough from the surface of the leptonic component, from which the particle is emitted, 
the particle is initially accelerated outward enough to escape from the Q-ball eventually. 

We show the solutions with various baryon numbers in Fig. 7. 2 The black circles indicates 
the solutions with B = L, which are realized at the Q-ball formation after the Affleck-Dine 
mechanism 3 . From there, the energy decreases along with L, which means that free particles 
are emitted and the Q-ball becomes electrically charged. The emission continues until the 
leptonic component vanishes for B <3x 10 4 . However, for B > 3 x 10 4 , the energy starts 
to increase in the middle of the decay, which means that the particle starts to form a cloud 
of bounded leptonic particles. In Fig. 8, we show the electric charge of Q-ball at which 
(dE/dL) B = 0. From the figure, ( dE/dL) B = 0 occurs for 

Q ~ 1.6 x 10 4 (71) 

2 The gaps in some data are due to the switching of algorithm for the computation. 

3 In the previous studies, the solutions were assumed to remain to be B = L, except for Ref. [14]. 
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Figure 5: 2-scalar gauged Q-balls in gauge mediation model with e 2 = 0.002, = 1 and 

B = 1.7 x 10 4 . 


and it is independent of baryon charge. This can be understood in the following way. When 
we approximate the energy of an emitted leptonic particle ( dE/dL)s as the energy with 
electricity switched off plus Coulomb energy, 


(dE 

\dL 


B 


~ CUql - 


- OJQL ~ 


e 2 Q 

AttRl 

e 2 Q 


(72) 


where we used the analytic expression with R = tt/co 0 . Equating ( dE/dL) B to zero, we 
obtain 

47T 2 

Q = — ~ 2.0 x 10 4 . (73) 

e 2 

This roughly explains Fig. 8. 

As in the case of 1-scalar gauged Q-ball, it is possible that 


(i dE/dB) L > 


(74) 
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Figure 6: 2-scalar gauged Q-balls in gauge mediation model with e 2 = 0.002, = 1 and 

B = 8.4 x 10 4 . 


due to the electrical repulsion. In order to check this, we only need to examine cui of each 
solution since ( dE/dB) L = uq as proven above. We illustrate where ( dE/dB) L = by 
black squares in Fig. 7. We plot the electric charge Q at which ( dE/dB) L becomes m,j, for 
each B in Fig. 9. The results can be fitted as 

log 10 Q = 3.0 + 3.0x lO-MogioS (75) 

This can be explained as follows. In the same way as deriving Eq. (72), the condition 


15 




























250000 


200000 


150000 


100000 


50000 


'v 


Q> 



®D, 


Ooo 




L 


®«(o] o ® 

) 0-® 


50000 


oo o © 

© 


d® 


100000 150000 

L 


O B=1.5xl0 5 
O B=1.0xl0 5 
O B=8.4xl0 4 
O B=5.6xl0 4 
O B=4.0xl0 4 
O B=3.0xl0 4 
O B=1.7xl0 4 
O B=L 

I I dE/dB=111(1, 
200000 250000 


Figure 7: Solutions with various baryon numbers in the E — L plane. The black circles are 
the solutions with B = L and the black squares are the solutions with (dE/8 B)l = m$. 


(' dE/dB )l = is written as 


Thus, we obtain 


= 


dE 

dB 


— UoB + 


— W0B + 


e 2 Q 
AttRb 
e%[ 

47r 2 


<^o B 


e 2 Q 


*^ U ° B 
_ i/4 


2y/2- 
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Q = EEb± b v 4 , 


V 0 1,4 e 2 


(76) 


(77) 


(78) 
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Figure 8: Electric charge Q when ( dE/dVjs = 0. 


which leads to 


logio Q = l°Sio 


f 2\/27rm$\ 

l +V ) 


+ I logio B 


~ 3.4+ 2.5 x 10 _1 log 10 S, 


(79) 


where we used the same approximation as before. This roughly explains our numerical 
results. 

The decay into protons is also expected to occur due to the electrical repulsion, even 
if we assume that the initial neutral Q-ball is stable against it. If this happens, baryonic 
component also decays off, and therefore charged Q-balls cannot be left in the universe. 
However, if the lcptonic decay stops before the electric charge becomes large enough for the 
baryonic decay to occur, the evolution stops and charged Q-balls may survive. One way to 
stop the lcptonic decay is that the leptonic cloud is close enough to the surface of the Q-ball 
so that the particle cannot come out. We derive the condition that this happens before the 
baryonic decay occurs. Let us suppose that the emitted leptons are electrons and roughly 
estimate the size of the leptonic cloud as Bohr radius. Then the condition that the cloud 


17 
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logB 


Figure 9: Electric charge when ( dE/dB)L = m$. 


radius becomes equal to the size of the Q-ball is written as 


47T 


m e Qe 2 


Q = 


= R= i=^ v " ll,Bl/a ’ 


m P e z 


This must happen for 


Q < 




(80) 

(81) 

(82) 


that is, before the baryonic decay starts to occur, where we have used the analysis parallel to 
that below Eq. (76). Therefore, the condition that the leptonic decay stops before baryonic 
decay occurs is given by 

V 0 


B > ~ 10 


30 


m 2 e m p 


(83) 


(10 6 GeV) 4 , 

This implies that if the baryon number is large enough, the leptonic decay may stop before 
the baryonic decay takes place, so that the charged Q-ball may survive as a relic in the 
universe. 

We see that for 

Vb 




7 T 4 m 4 


10 


35 


(10 6 GeV) 4 


(84) 
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Bohr radius is already smaller than the Q-ball size when the cloud is about to be formed, 
which means that the evolution stops without forming the cloud. Thus, more accurately, if 


m 2 m 2 7r 4 mf ’ 



(85) 


then the charged Q-balls are expected to survive with a cloud surrounding it, and if 



( 86 ) 


then the charged Q-balls are expected to survive without the cloud 1 . 


Finally, the dashed line in Fig. 7 illustrates Q m ax f° r each B, above which Q-ball solutions 
cannot exist. This results from electrical repulsion due to the gauge field. 

5 Conclusions and discussion 

In this paper, we considered gauged Q-balls in the two scalar model in order to discuss 
the evolution of neutral Q-balls which are formed from the flat direction in the Affleck- 
Dine mechanism and the possibility of realization of gauged Q-balls during their evolution. 
We approximated the evolution as a sequence of charged Q-ball solutions, which is implied 
from the situation that only the leptonic component decays off. As a result, a Coulomb 
potential arises and the Q-ball becomes electrically charged as expected. In other words, it 
is energetically favored for leptonic decay to occur. However, since there is an upper bound 
on charge of gauged Q-ball, the amount of decay is limited as well, which we examined 
quantitatively. In addition, if the baryon number of the initially formed Q-ball is large 
enough, the electric charge of the Q-ball grows enough so that the particle which is emitted 
from the Q-ball is bound to it. The baryonic decay is also expected to occur by virtue of 
the electrical repulsion, which leads to the vanishing of the charged Q-balls in the universe. 
However, it is expected that if the leptonic cloud is close enough to the surface of the Q-ball, 
the leptonic decay can stop before the electric charge becomes large enough for the baryonic 
decay to occur, so that the evolution stops and charged Q-balls can survive. We roughly 
estimated when this can happen, and as a consequence, we found that there exists a lower 
bound on the baryon number. 

Suppose that dark matter consists of Q-balls. Since a Q-ball is known to absorb protons 
and emit pions, the neutral Q-balls can be detected by Super-Kamiokande [9, 10, 11] or 

4 If B > 4Vg/i r 4 ?n 4 , the baryonic decay starts to occur when Q = 4m p /m e e 2 , and this is larger than 
Q = 47r 2 /e 2 , which is when the evolution is expected to stop. Thus, there is no need to worry that the 
baryonic decay occurs before the evolution stops. 
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IceCube [19], while these detectors are not suited for detection of the charged Q-balls since 
charged Q-balls cannot absorb the protons due to the electrical repulsion. Thus, the electric 
charge is expected to make differences in the experimental signatures of the relic Q-balls. 
The charged Q-balls are expected to behave as some kind of nuclei, and are known to be 
detectable by such detectors as MACRO [9, 12] and the observational bounds on mass and 
flux of the relic charged Q-balls are obtained [10]. The emitted particles are also expected 
to contribute to the energy components of the universe and lcptonic particles must satisfy 
the observational bounds on neutrino component from WMAP [20], which must be verified 
in future work. 

Although we assumed that the decay occurs if energetically allowed, we did not specify 
which kind of decay we are considering, for instance species of particle which is emitted, an 
actual decay rate etc. These informations must be added when we consider electrodynamical 
effects and actually solve time development of Q-ball. Lastly, here we are considering the 
simplest model of a flat direction, but in reality, we may need to consider more complex flat 
directions, possibly in nonabclian gauge theories as well. 
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